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Abstract. We report the discovery of WASP-6b, an inflated sub- Jupiter mass planet transiting every 
3. 36100601q 0000035 days a mildly metal-poor solar-type star of magnitude V=11.9. A combined analysis of 
the WASP photometry, high-precision followup transit photometry and radial velocities yield a planetary mass 
M p = 0.503±g g|| Mj and radius R v = 1.224^;^ Rj, resulting in a density p p = 0.27 ± 0.05 pj. The mass 
and radius for the host star are M* = 0.88l^o8 

M Q and R, = 0.8701^.^36 Rq. The non-zero orbital eccentricity 
e = 0.0541q'qj 8 that we measure suggests that the planet underwent a massive tidal heating ~ 1 Gyr ago that 
could have contributed to its inflated radius. High-precision radial velocities obtained during a transit allow us 
to measure a sky-projected angle between the stellar spin and orbital axis (5 — llllg deg. In addition to similar 
published measurements, this result favors a dominant migration mechanism based on tidal interactions with a 
protoplanetary disk. 

Key words, binaries: eclipsing - stars: individual: WASP-6 - planetary systems - techniques: photometric - 
techniques: radial velocities - techniques: spectroscopic 



1. Introduction 

Transiting planets play an important role in our under- 
standing of the nature of the extrasolar planetary objects. 
They are the only exoplanets for which an accurate mea- 
surement of the mass and radius is available. The deduced 
density is a key parameter to constraint theoretical mod- 
els for the formation, evolution and structure of planets 
(e.g. Fortney et al. 2007; Liu et al. 2008). For the bright- 
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est transiting systems, a study of the atmospheric compo- 
sition and physics is possible, even with existing instru- 
ments like HST or Spitzer (e.g. Charbonneau et al. 2008; 
Swain et al. 2008). The discovery rate of transiting plan- 
ets has increased recently thanks mainly to the efficiency 
of the CoRoT space-based survey (Baglin et al. 2006) and 
of a handful of ground-based wide-field surveys target- 
ing rather bright stars (V < 13): HATNet (Bakos et al. 
2004), WASP (Pollaco et al. 2006), TrES (O 'Donovan et 
al. 2006), and XO (McCullough et al. 2005). 

The ~ 50 transiting planets known at the time of this 
writing show a broad range of mass and radius. Their 
masses go from 23 Af® for the hot Neptune GJ436b 



2 



M. Gillon et al.: The transiting sub- Jupiter mass planet WASP-6b 



(Butler et al. 2004; Gillon et al. 2007) to more than 10 Mj 
for XO-3 (Johns-Krull et al. 2008). Many planets have a 
size in concord with basic models of irradiated planets (e.g. 
Burrows et al. 2007, Fortney et al. 2007), some of them 
like HD 149026 b (Sato et al. 2005) appearing to be very 
rich in heavy elements. Nevertheless, a few planets like 
HD 209458b (Charbonneau et al. 2000, Henry et al. 2000) 
are 'anomalously' large. Several hypothesis have been pro- 
posed to explain this radius anomaly, most importantly 
tides (Bodenheimer et al. 2001; Jackson et al. 2008b), tides 
with atmospheric circulation (Guillot & Showman 2002) 
and enhanced opacities (Guillot et al. 2006, Burrows et 
al. 2007). The existence of several correlations between 
parameters of transiting systems has been proposed, for 
instance between the planet mass and the orbital period 
(Mazeh et al. 2005; Gaudi et al. 2005) and between the 
heavy-element content of the planet and the stellar metal- 
licity (Guillot et al. 2006; Burrows et al. 2007). The astro- 
physics supporting these correlations has still to be fully 
understood. 



It is highly desirable to detect and characterize thor- 
oughly many more bright short period transiting systems 
to improve our understanding of the highly irradiated 
gaseous planets and to constraint the structure and evo- 
lution models for these objects. With its very high detec- 
tion efficiency, the WASP transit survey is making a large 
contribution to this goal. It is the only transit survey op- 
erating in both hemispheres: it uses an instrument named 
WASP-North and located at La Palma to search for plan- 
ets from the Northern hemisphere and a twin instrument 
named WASP-South and located at Sutherland to do the 
same from the Southern hemisphere. Each of these instru- 
ments covers a huge field of view of 482 square degrees per 
pointing, allowing them to search for transiting planets in 
a large portion of the sky. Due to the brightness of the 
host stars, planets detected by WASP are very good tar- 
gets for high-precision followup observations. For instance, 
it is possible to measure for most of them the alignment 
between the stellar rotation axis and the planetary orbital 
axis via the observation of the Rossiter-McLaughlin effect 
(RM; Queloz et al. 2000). The measured value for this 
spin-orbit angle is a strong constraint for inward planetary 
migration models (see Winn 2008 and references therein) . 

We report here the discovery and characterization of 
WASP-6b, a new sub-Jupiter mass planet transiting a 
mildly metal-poor solar- type star of magnitude V=11.9. 
We present in Section 2 the WASP discovery photometry 
plus high precision followup transit photometry and radial 
velocity measurements confirming the planetary nature of 
WASP-6b and including the observation of a spectroscopic 
transit. Section 3 presents the determination of the host 
star parameters. Our determination of the system param- 
eters is presented in Section 4. These parameters are dis- 
cussed in Section 5. 




Phase 



Fig. 1. WASP photometry of WASP-6 phase-folded with 
the best-fit period from the transit search algorithm pre- 
sented in Collier Cameron et al. (2006). 

2. Observations 

2.1. WASP photometry 

The host star 1SWASP J231237.75-224026.1 (= USNO- 
B1.0 0673-1077008 = 2MASS 23123773-2240261; here- 
after WASP-6) was observed by WASP-South during 
the 2006 and 2007 observing seasons, covering the in- 
tervals 2006 May 07 to 2006 November 12 and 2007 
July 05 to 2007 November 13 respectively. The 9630 
pipeline-processed photometric measurements were de- 
trended and searched for transits using the methods de- 
scribed in Collier Cameron et al. (2006). The selection 
process (Collier Cameron et al. 2007) elected WASP-6 as 
a high priority candidate presenting a periodic transit-like 
signature with a period of 3.361 days. A total of 18 tran- 
sits are observed in the data. Figure 1 presents the WASP 
photometry folded with the best-fit period. 

2.2. High-S/N transit photometry 

Followup transit photometry was obtained on 2007 
October 13 using the 2048x2048 pixel 2 camera 
HawkCam2 (Wilson et al. 2008, Anderson et al. 
2008) on the 2.0-m Faulkes Telescope South (FTS) at 
Siding Spring Observatory . The camera has a scale of 
0.135 arcseconds/pixel and a field of view of ~ 4.6 x 4.6 
arcminutes 2 . We observed the target field using the SDSS 
V band in the 2x2 bin mode to improve the duty cycle. 
We acquired 247 frames of 60 sec exposure during the 
run. The telescope was sufficiently defocussed to keep 
the stellar flux within the linear range of the CCD. 
The images were bias subtracted and flat-field corrected 
with a master bias and twilight flat field images using 
IRAF0. DA0PH0T aperture photometry (Stetson 1987) was 
performed around the target and comparison stars. We 
substracted a linear fit from the differential magnitudes 
as a function of airmass to correct for the different 
colour dependance of the extinction for the target and 

1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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Photometry vs Phase 



1.01 
1.00 

0.99 

u_ 

0.98 
0.97 
0,96 

-0.02 

-0.01 

= 0.00 
Li. 

0.01 

0.02 




000 

Phase 
Photometry vs Phase 



1.01 
1.00 

0,99 

i 

0.9S 





Fig. 2. FTS i'-band (top) and LT/RISE V + R (bottom) 
transit photometry for WASP-6 and residuals after sub- 
traction of the best-fit transit curve (superimposed in 
blue). 

comparison stars. The linear fit was calculated from the 
out-of-transit (OOT) data and applied to all the data. 
The corresponding fluxes were then normalized using 
the OOT part of the photometry. We discarded the first 
17 measurements because they were obtained during 
twilight. Fig. 2 shows the resulting lightcurve folded on 
the best-fit orbital period and the residuals obtained 
after removing the best-fit transit model (see Section 
4). Their rms is 1.67 x 10~ 3 . This can be compared to 
9.54 x 1(T 4 , the mean theoretical error bar taking into 
account photon, read-out, scintillation and background 
noises. 

High precision transit observations of WASP-6 were 
also carried out using the 1024x 1024 pixel 2 thcrmoelectri- 
cally cooled frame transfer CCD camera RISE mounted on 
the 2-m Liverpool Telescope (LT) in La Palma (Steele et 
al. 2008). The camera has a scale of 0.55 arcseconds/pixel 
and a total field of view of ~ 9.4 x 9.4 arcminutes 2 . We 



observed the target field using a single broad band V + R 
filter in the 2x2 bin mode. We acquired 4200 frames of 3 
sec exposure on the night of 2008 July 25 and 2880 frames 
of 5 sec exposure on the night of 2008 August 11. The 
telescope was, here too, defocussed. A similar reduction 
procedure as for the FTS photometry was used. The re- 
sulting normalized light curves of WASP-6 folded with the 
best- fit orbital period are shown in Fig. 2. The rms of the 
residuals is respectively 0.54 % and 0.5 % for the first and 
second run, while their mean theoretical error bar are 0.51 
% and 0.40 %. 



2.3. Spectroscopy 

As soon as WASP-6 was identified as a high priority tar- 
get, spectroscopic measurements were obtained using the 
CORALIE spectrograph mounted on the Euler Swiss tele- 
scope (La Silla, Chile) to confirm the planetary nature of 
the eclipsing body and measure its mass. WASP-6 was ob- 
served from 2007 September 16 to 2007 October 26 and 
from 2008 September 11 to 2008 September 25. Radial ve- 
locities (RV) were computed by weighted cross-correlation 
(Baranne et al. 1996; Pepe et al. 2005) with a numerical 
G2-spectral template. RV variations of semi-amplitude ~ 
75 m s _1 were detected consistent with a planetary-mass 
companion whose period closely matches that from the 
WASP transit detections. 

44 additional spectroscopic measurements were ob- 
tained with the HARPS spectrograph (Mayor et al. 2003) 
based on the 3.6-m ESO telescope (La Silla, Chile) in 
the context of the programs 082.C-0040(E) and 082. C- 
0608(E). These programs aim to improve the characteriza- 
tion of WASP transiting planets. As CORALIE, HARPS 
is a cross-dispersed, fiber-fed, echelle spectrograph dedi- 
cated to high-precision Doppler measurements. HARPS 
data were reduced with a pipeline very similar to the 
CORALIE one. In addition to several measurements cov- 
ering the whole orbital phase, high-cadence measurements 
of a spectroscopic transit were obtained with HARPS on 
2008 October 08 in order to determine the sky-projected 
angle between the planetary orbital axis and the stellar 
rotation axis and included two points taken the night be- 
fore, a point as far as possible from the transit on the 
transit night and a point the night after. This strategy 
aims to determine the systematic RV with greater accu- 
racy than if the RM effect was taken on its own, assuming 
that stellar activity is the same over the three nights. 

Our RV measurements are listed in Table 1 
(CORALIE) and Table 2 (HARPS) and arc shown phase- 
folded and over-plotted with the best-fitting orbital+RM 
model in Fig. 3. 

To exclude that the RV signal shown in Fig. 3 is due 
to spectral line distortions cause by a blended eclipsing 
binary, the CORALIE and HARPS cross-correlation func- 
tions were analyzed using the line-bisector technique de- 
scribed in Queloz et al. (2001). No evidence for a corre- 
lation between the bisector spans and the RV variations 
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Fig. 3. Top: The RV measurements of WASP-6 ob- 
tained with CORALIE (red triangles) and HARPS (green 
squares). The systematic velocity has been subtracted. 
The solid line is the MCMC solution (see Section 4) ; it in- 
cludes the RM effect. Bottom: zoom on the transit phase 
showing the RM effect. 

was found (Fig. 4). The most likely cause for the periodic 
signal observed in photometry and RV measurements and 
for the RM effect observed on 2008 October 08 is thus 
the presence of a giant planet transiting the star WASP-6 
every 3.36 days. 

3. WASP-6 Stellar Parameters 

The individual CORALIE and HARPS spectra are rela- 
tively low signal-to-noise, but when co-added into O.OlA 
steps they give a S/N of in excess of 100:1 which is suit- 
able for a photospheric analysis of WASP-6. The standard 
pipeline reduction products were used in the analysis. 

The analysis was performed using the uclsyn spectral 
synthesis package (Smith 1992; Smalley et al. 2001) and 
ATLAS9 models without convective overshooting (Castelli, 
Gratton & Kurucz 1997). The H a line were used to de- 
termine the effective temperature (T c g), while the Na I D 
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Fig. 4. Bisector versus RV measured from all the observed 
CORALIE (top) and HARPS (bottom) spectra. We adopt 
uncertainties of twice the RV uncertainty for all bisector 
measurements. There is no correlation between these two 
parameters indicating the RV variations are not caused by 
stellar activity or line-of-sight binarity. 



and Mg I b lines were used as surface gravity (log g) di- 
agnostics. The parameters obtained from the analysis are 
listed in Table E] 

The equivalent widths of several clean and unblended 
lines were measured. Atomic line data was mainly taken 
from the Kurucz & Bell (1995) compilation, but with up- 
dated van der Waals broadening coefficients for lines in 
Barklem et al. (2000) and log gf values from Gonzalez & 
Laws (2000), Gonzalez et al. (2001) or Santos et al. (2004). 
A value for microturbulence (£t ) was determined from Fe I 
using Magain's (1984) method. The ionization balance be- 
tween Fe I and Fe n and the null-dependence of abundance 
on excitation potential were used as an additional the T e g 
and logg diagnostics (Smalley 2005). 

We have determined the elemental abundances of sev- 
eral elements (listed in Table [3]) from their measured 
equivalent widths. The quoted error estimates include that 
given by the uncertainties in T e g, logg and £ t > as well as 
the scatter due to measurement and atomic data uncer- 
tainties. In our spectra the Li I 6708 A line is not detected 
(EW < 2mA) , allowing us to derive an upper-limit on the 
Lithium abundance of log n(Li/H) + 12 < 0.5. The lack 
of lithium implies an age in excess of ~3 Gyr (Sestito & 
Randich 2005). 

Projected stellar rotation velocity (V ro tsmI) was de- 
termined by fitting the profiles of several unblended Fe I 
lines in the HARPS spectra. We used a value for macro- 
turbulence (w mac , see Gray 2008) of 2 km s _1 and an in- 
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BJD-2,400,000 RV a RV BS 



(days) (km s 1 ) (km s 1 ) (km s 1 ) 



54359. 


687716 


11 


.48815 





.01636 


0.01618 


54362. 


.573582 


11 


.43618 





.01602 


-0.05289 


54364, 


628975 


11 


.44037 





.02006 


-0.02655 


54365. 


707645 


11 


.42032 





.01334 


-0.03610 


54372. 


735312 


11 


.42022 





.01288 


-0.06013 


54377. 


716124 


11 


.48411 





.01452 


-0.04330 


54377. 


.739618 


11 


.48861 





.01393 


-0.04916 


54378, 


693586 


11 


.38638 





.01271 


-0.03449 


54378, 


719349 


11 


.38469 





.01306 


0.01096 


54379, 


690005 


11 


.45960 





.01414 


-0.03955 


54379, 


713452 


11 


.46897 





.01333 


-0.04065 


54380, 


563292 


11 


.53322 





.01202 


-0.01290 


54380, 


586751 


11 


.51439 





.01142 


-0.02735 


54382, 


.617808 


11 


.39778 





.01422 


-0.00658 


54383, 


601838 


11 


.50697 





.01225 


-0.03398 


54383, 


625413 


11 


.49114 





.01253 


-0.04571 


54385, 


740623 


11 


.36965 





.02555 


-0.10750 


54386, 


638647 


11 


.46303 





.01407 


-0.01829 


54386, 


.664860 


11 


.49377 





.01384 


-0.02509 


54387, 


668883 


11 


.48742 





.02578 


-0.06456 


54387, 


690688 


11 


.49402 





.03290 


0.04418 


54390, 


521372 


11 


.52130 





.01273 


-0.08028 


54390, 


544901 


11 


.51282 





.01314 


-0.03094 


54398, 


.647981 


11 


.36002 





.01361 


-0.02246 


54398, 


671428 


11 


.38428 





.01373 


-0.02346 


54399, 


.616495 


11 


.41113 





.01204 


-0.05097 


54399, 


652117 


11 


.42110 





.01225 


-0.05205 


54720, 


614742 


11 


.50358 





.01721 


-0.05665 


54722, 


.651535 


11 


.45913 





.01483 


-0.04094 


54724, 


579899 


11 


.40131 





.01616 


-0.00673 


54725, 


545387 


11 


.43430 





.01835 


-0.06935 


54726, 


613463 


11 


.52363 





.01670 


-0.01711 


54730, 


622609 


11 


.47459 





.01459 


-0.00302 


54732, 


703073 


11 


.47021 





.01808 


-0.03689 


54734, 


747409 


11 


.41662 





.01649 


0.00588 



Table 1. CORALIE radial velocity measurements for 
WASP-6 (BS = bisector spans). 



strumental FWHM of 0.060 ± 0.005 A, determined from 
the telluric lines around 6300A. A best fitting value of 
V^otsin/ = 1.4 ± 1.0 km s _1 was obtained. If, however, 
mactroturbulence is lower, then higher rotation values are 
found, with V^. ot sinJ = 3.0 ± 0.5 km s _1 obtained for 
w mac = km s . If, on the other hand, w mac is slightly 
higher than 2 km s _1 , then it is possible that V rot sin I is 
close to, or even, zero. 

In addition to the spectral analysis, we have also used 
broad-band photometry from TYCHO-2, USNO-B1.0 R- 
mag., CMC14 r', DENIS and 2MASS to estimate the to- 
tal observed bolometric flux. The Infrared Flux Method 
(Blackwell & Shallis 1977) was then used with 2MASS 
magnitudes to determine T e g and stellar angular diameter 
(9). This gives Tcff^ 5470 ± 130 K, which is in close agree- 
ment with that obtained from the spectroscopic analysis 
and implies a spectral type of G8V (Gray 2008). 



BJD-2,400,000 ~ RV a RV BS 



(days) (km s 1 ) (km s x ) (km s x ) 



54746. 


564605 


11. 


.54738 


0. 


,00576 


0.00133 


54746. 


750752 


11. 


.54840 





,00555 


-0.01787 


54747. 


504842 


11. 


.51620 





00559 


-0.02660 


54747. 


579780 


11 


.50597 





,00711 


-0.03234 


54747. 


586192 


11 


.51299 





,00677 


0.00139 


54747. 


592603 


11. 


.49959 





,00665 


-0.00608 


54747. 


598784 


11. 


.50313 





,00781 


-0.01701 


54747. 


605195 


11 


.50141 





,00819 


-0.07376 


54747. 


611549 


11. 


.48940 





,00885 


-0.06980 


54747. 


617718 


11 


.49653 





,00811 


-0.03563 


54747. 


624372 


11. 


.48137 





,00808 


-0.00962 


54747. 


630668 


11. 


.50486 





,00716 


-0.02661 


54747. 


636837 


11. 


.51929 





,00812 


0.02936 


54747. 


643307 


11 


.49737 





,00808 


-0.01246 


54747. 


650089 


11. 


.52399 





,00766 


-0.00152 


54747. 


655956 


11. 


.50607 





,00612 


-0.01990 


54747. 


662310 


11. 


.50622 





,00648 


-0.00217 


54747. 


668664 


11. 


.49528 





,00706 


-0.05534 


54747. 


675076 


11. 


.48567 





,00634 


-0.01449 


54747. 


681372 


11. 


.49207 





,00624 


0.02179 


54747. 


687841 


11. 


.46478 





,00616 


-0.03649 


54747. 


694137 


11 


.45835 





,00571 


-0.05455 


54747. 


700491 


11 


.47240 





,00618 


-0.01938 


54747. 


706729 


11. 


.46878 





,00655 


-0.00382 


54747. 


713256 


11. 


.45266 





,00619 


0.02667 


54747. 


719494 


11. 


.46519 





00625 


-0.03618 


54747. 


725906 


11 


.46526 





,00591 


0.00022 


54747. 


732260 


11. 


.48427 





,00603 


0.00930 


54747. 


738672 


11. 


.46037 





,00689 


0.01768 


54747. 


744968 


11. 


.46863 


0, 


,00668 


0.00392 


54747. 


751264 


11. 


.45986 


0, 


,00685 


-0.02154 


54747. 


.757617 


11. 


.46432 





,00734 


0.01274 


54747. 


764029 


11. 


.47662 





,00759 


-0.04832 


54747. 


770383 


11 


.46143 





,00783 


-0.07680 


54747. 


776621 


11. 


.49690 





,00818 


0.03016 


54747. 


783206 


11. 


.46028 





,00886 


-0.03313 


54747. 


789329 


11 


.47787 





,00875 


-0.02263 


54747. 


795798 


11. 


.45389 





,00916 


0.01502 


54748. 


775349 


11. 


.42689 





,00406 


-0.02922 


54749. 


703663 


11. 


.53125 





,00270 


-0.02405 


54750. 


655365 


11. 


.53933 





,00240 


-0.02941 


54753. 


741910 


11. 


.54705 





,00414 


-0.01725 


54754. 


750340 


11. 


.43295 





,00299 


-0.02203 


54755. 


730495 


11. 


.43908 





,00261 


-0.02337 



Table 2. HARPS radial velocity measurements for 
WASP-6 (BS = bisector spans). 



4. Derivation of the system parameters 

We derived stellar and planetary parameters for the 
system by fitting simultaneously the WASP, FTS and 
LT/RISE photometry with the CORALIE and HARPS 
RVs. These data were used as input into the Markov 
Chain Monte Carlo (MCMC; Ford 2006) code described 
in Gillon et al. (2008). MCMC is a Bayesian inference 
method based on stochastic simulations and provides the 
a posteriori probability distribution of adjusted param- 
eters for a given model. Here the model is based on a 
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Parameter 



Value 



R.A. (J2000) 
Dec (J2000) 
V 

Toff- 
log 9 

ft 

Vrot sin I 



23 ft 12 m 37.74" 
-22°40'26".2 
11.9 

5450 ± 100 K 
4.6 ± 0.2 
1.0 ± 0.2 km s" 1 
1.4 ± 1.0 km s _1 



[Na/H] 


-0. 


17 


± 


0.06 


[Mg/H] 


-0. 


.13 


± 


0.07 


[Al/H] 


-0. 


15 


± 


0.10 


[Si/H] 


-0. 


12 


± 


0.08 


[Ca/H] 


-0. 


09 


± 


0.10 


[Sc/H] 


-0. 


22 


± 


0.15 


[Ti/H] 


-0. 


05 


± 


0.09 


[V/H] 


-0. 


02 


± 


0.08 


[Cr/H] 


-0. 


17 


± 


0.09 


[Mn/H] 


-0. 


20 


± 


0.13 


[Fe/H] 


-0. 


20 


± 


0.09 


[Co/H] 


-0. 


16 


± 


0.14 


[Ni/H] 


-0. 


21 


± 


0.08 



log N(Li) 
T eff (IRFM) 
0(IRFM) 



< 0.5 
5470 ± 130 K 
0.037 ± 0.002 mas 



Table 3. Stellar parameters for WASP-6. 



star and a transiting planet on a keplerian orbit about 
their center of mass. Specifically, we used the photometric 
transit model of Mandel & Agol (2002) and the spectro- 
scopic transit model of Gimenez (2006) in addition to a 
classical Keplerian model for the orbital part of the RV 
variations. To model the transit lightcurves, a quadratic 
limb darkening law was assumed, with coefficients interpo- 
lated from Claret's tables (2000; 2004) for the appropriate 
photometric filters. For the RISE broad band filter, the av- 
erage from V and R bands was taken to be our theoretical 
limb darkening parameters. 

We used 16 jump parameters in our MCMC simula- 
tions: the orbital period P, the time of minimum light 
To, the transit depth D, the total transit width W, the 
impact parameter b, the stellar mass M», the orbital RV 
semi-amplitude K, a systematic radial velocity 7 for each 
spectroscopic instrument (HARPS and CORALIE), the 
parameters ecosw and esinw where e is the orbital ec- 
centricity and lo is the argument of periastron, the prod- 
ucts V rot sin / cos (3 and V rot sin I sin j3 where V rot sin / is 
the projected stellar rotational velocity and f3 is the spin- 
orbit angle (see Gimenez 2006), and a normalization factor 
for each of the 4 light curves (assuming the same normal- 
ization for the whole SW photometry). As explained in 
the now abundant literature on the application of MCMC 
to perform Bayesian inference for transiting planets (see 
e.g. Collier Cameron et al. 2007 and references therein), 
each MCMC simulation is composed of a large number of 



consecutive steps for which the jump parameters are ran- 
domly modified or not depending of the result of a test on 
the merit function (MF). The MF used here is the sum 
of the x 2 f° r au the data with respect to the models added 
to a Bayesian prior on V rot sin I and M* representing our 
constraints on these parameters from spectroscopy: 



MF = x 2 - 



(V rot sin I - (V r 



a 



V rot sin / 



^tsinPjo) 2 | (M„ - (M«)q) 



where (V^ tsin/)o — 1.4 km s 1 , <rv rot sini = 1 km 
s _1 , M* — 0.87 and um, — 0.08. These last two values 
were obtained by interpolation of Girardi stellar evolution 
models (Girardi et al. 2000) in order to find the mass 
and age that best match the spectroscopic parameters. 
We notice that our data do not constrain strongly M* and 
that it is a free parameter under the control of a Bayesian 
prior in our simulations only to propagate its uncertainty 
to the other physical parameters. 

A first MCMC run was performed and led to a refined 

1/3 

value for the stellar density. We converted it to R^/Mj 
in solar units, and compared this property and the stel- 
lar temperature to the Girardi models interpolated at -0.2 

1 /3 

metallicity. The quantity, R^jMj , depends only on the 
observed transit properties (duration, depth, impact pa- 
rameter, and orbital period) and is independent of the 
measured temperature. We generated the same property 
from the mass and log g values in the models, and then 
interpolated the models in the R/M 1 / 3 -^^ plane to deter- 
mine a mass and age for WASP-6. We interpolated linearly 
along two consecutive mass tracks to generate an equal 
number of age points between the zero-age main sequence 
and the evolutionary state where the star reaches the end 
of core hydrogen burning. We then interpolated between 
the mass tracks along equivalent evolutionary points to 
find the mass and age from the models that best match 
the stellar density derived from the MCMC and the effec- 
tive temperature. In this way, we obtained a value for the 
stellar mass of, M* = 0.83io.o9 ^0 ari( l a derived age for 
the system of 11 ± 7 Gyr. 

The best-fitting model found in the first MCMC run 
was used to estimate the level of correlated noise in each 
photometric time-series and a jitter noise in the RV time 
series. For each photometric time-series, the red noise was 
estimated as described in Gillon et al. (2006), by com- 
paring the rms of the unbinned and binned residuals. We 
used a bin size corresponding to a duration of 25 minutes, 
similar to the timescale of the ingress/egress of the tran- 
sit. For the SW data, the red noise was estimated to be 
negligible when compared to the theoretical error bar of 
the measurements and it was thus neglected. The deduced 
red noise values (Table 4) were added quadratically to 
the theoretical uncertainties of each corresponding time- 
series. No jitter is detected in the CORALIE data. For the 
HARPS data, a significant jitter is obtained, but it seems 
to be originating mostly from the residuals of the RM ef- 
fect and is probably more due to lower-than-usual S/N 
on the spectra and a worsening of airmass (reaching 1.8 
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Fig. 5. R/M 1 ' 3 (in solar units) versus effective tempera- 
ture for WASP-6 compared to the theoretical stellar stel- 
lar evolutionary models of Girardi et al. (2000). The la- 
beled mass tracks are for 0.8, 0.9 and 1.0 Mq and the 
isochrones are 100 Myr (solid), 1 Gyr (dotted), 5 Gyr (dot- 
dashed), 10 Gyr (dotted). We have interpolated the tracks 
at -0.2 metallicity and have included the uncertainty on 
the metallicity (+- 0.1) in the overall uncertainties on the 
mass and the age. According to the models, the host star 
has an age of 11 ± 7 Gyr. 

at the end of the sequence) than stellar activity. For this 
reason, no jitter noise was added to the RV uncertainties. 

Using the updated value of the stellar mass as initial 
value, a second MCMC run was then performed. This 
chain allowed a large safety burn-in period of discarded 
50,000 steps followed by a simulation of 500,000 steps al- 
lowing a robust determination of the a posteriori prob- 
ability distributions for the jumped parameters. The pa- 
rameters set (jump + deduced parameters) corresponding 
to the lowest MF was considered as the best solution, and 
for each parameter upper and lower 1-rj error bars were 
obtained from respectively the 68.3% larger and smaller 
values. Best-fitting jump + physical parameters are shown 
in Table 5. The reduced x 2 of the best-fitting solution is 
0.86. 

5. Discussion 

The large radius of WASP-6b (~ l.2Rj) and the metal de- 
ficiency of its host star strengthen the existence of a corre- 
lation between the heavy-element content of giant planets 
and the stellar metallicity (Guillot et al. 2006; Burrows 
et al. 2007). With half of the mass of Jupiter and a ra- 
dius significantly larger, WASP-6b appears nevertheless 
too large for basic models of irradiated planets (Burrows 
et al. 2007a; Fortney et al. 2007), even if an absence of 
core is assumed. For instance, tables presented in Fortney 
et al. (2007) predict a maximum radius of ~ l.lRj for a 
0.5 Jupiter-mass planet orbiting at 0.045 AU of a 4.5 Gyr 
solar-type star. WASP-6 is smaller, cooler and probably 



older than the Sun, so \.2Rj is clearly too large for these 
models. In this context, it is worth noticing the non-null 
eccentricity that we infer for its orbit (e = 0.054+qq^). 
The fact that the planetary orbit is still not circularized 
despite the large age of the system indicates that the tidal 
evolution of WASP-6b probably played an important role 
in its energy budget. As outlined by Jackson et al. (2008b), 
tidal heating could have been large enough for many close- 
in planets to explain at least partially the large radius of 
some of them. To assess the past and future tidal evolu- 
tion of WASP-6b, we integrated the equations for da/dt 
and de/dt presented in Jackson et al. (2008a) and com- 
puted at each step the tidal heating rate H using the for- 
mula presented in Jackson et al. (2008b). We assumed 
values of Q' p = 10 6 5 and Q'„ = 10 5 5 for respectively the 
planetary and stellar tidal dissipation parameter^. These 
values were found by Jackson et al. (2008a) to concili- 
ate the eccentricity distribution of close-in planets before 
their tidal evolution to the one of the planets detected 
farer from their star. We also took into account the evo- 
lution of the stellar rotation period due to the tide raised 
by the planet using (Goldreich & Soter 1966): 

dQ 9 R 3 M 2 

"-'»<"• -4 G Jm^- < 2 > 

where G is the gravitational constant, n is the mean or- 
bital motion, is the stellar spin angular rate and a* = 
I*/(M*R 2 ) with /* being the moment of inertia though 
the spin axis of the star. For a», we assumed a value of 
0.07 (Patzold et al. 2004). To assess the reliability limits 
of the model, we also computed the evolution of the total 
angular momentum of the system (assuming a negligible 
contribution of the planet rotation): 

Uot = ^ UM ! lt na 2 Vl-e 2 + a^R 2 ^ , (3) 
M* + Mp 

Neglecting the possible decrease due stellar wind (Dobbs- 
Dixon et al. 2004), L to t should be a conserved quantity 
during the whole tidal evolution of the system. 

Fig. 6 shows the obtained evolution for a, e, H, L to t 
and the orbital and stellar rotation period from 2 Gyr ago 
to 5 Gy in the future. Interestingly, the model predicts 
(1) that the eccentricity and semi-major axis of WASP-6b 
were significantly larger in the past, (2) that the orbit will 
be fully circularized one Gy from now, and (3) that the 
planet will continue to slightly approach the star until fi- 
nally reaching its Roche limit. This last results agrees well 
with the fact that the ratio L tot /L c , where L c is critical 
angular momentum (see Levrard et al. 2009), has a value 
of ~ 0.6, implying that the system is tidally unstable and 
will ultimately merge. Levrard et al. (2009) showed that 
all the other transiting systems, except HAT-P-2, are in 
the same case. 

2 We use here the same convention than Jackson et al. 
(2008a): the coefficients Q' p and Q' t used here are equal to 
the actual tidal dissipation parameters Q p and Q* multiplied 
by the ratio 3 /2k where k is the Love number. 
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Photometric time-series 


Red noise [ppm] 




r 1 o 


040 




RlbE-1 


317 




RISE-2 


770 




RV time-series 


Jitter [ms _1 ] 


systematic RV [km s 1 ] 


CORALIE 





11.449 


HARPS 


6.4 


11.485 



Table 4. Deduced values for the photometric red noise (top), and the RV jitter and systematic velocities (bottom). 



Parameter 



Value 



Units 



Transit epoch To 

Orbital period P 

Planet/star area ratio (R p /R s ) 2 

Transit duration tr 

Impact parameter b 

RV semi-amplitude K 

ecoscj 

esincj 

Vrot sin / cos f3 
V r ot sin I sin (3 



2454596.43267to;°ooio 
3 36100 60+ 0000022 

O.OOlUUOU_ .0000035 

02092+ 00019 

U.UiU3i_ .00025 



0.1086018S 
0.261 
74.31 
-0.007! 



~+0.07 
-0.11 
3+1.7 
-1.4 
7 +0.()ll 
-0.008 

0.054t°;°l? 



+0.28 

J-- ' -0.10 

q on+0.49 
o.oz_ 50 



HJD 

days 

days 
R* 



Orbital semi-major axis a 
Orbital inclination i 
Orbital eccentricity e 
Argument of periastron u> 
Spin-orbit angle (3 



n491 +u uuus 

U.U4Zl_o .0013 

88.47l^| 

U.UO4_ oi 5 

i 7n+ - 12 

L - ' u -0.23 
9D+ - 25 



AU 
degrees 

rad 
rad 



Stellar mass M» 

Stellar radius i?* 

Stellar surface gravity log g, 

Stellar density p* 

Projected rotational velocity V ro tsmI 



0.880±g;g|g 

0.8701°;^ 
4.50 ±0.06 

i-6o±8:!? 



M 
J?0 



P© 
km s _1 



Planet radius i? p 

Planet mass M p 

Planetary surface gravity log g p 

Planet density p v 

Planet temperature (A = 0,f= 1/4) T eff 



1.224" 



0.052 



U.OUO_ 038 

7.857 ±0.028 
0.27 ±0.05 
119411? 



Rj 
M., 
[cgs] 

PJ 
K 



Table 5. WASP-6 system parameters and \-a error limits derived in this work. 



Under this tidal evolution model, WASP-6b was 
brought to a distance > 0.05 AU of its host star in the 
very early life of the system, then its orbital evolution has 
been totally dominated by tides until now. This evolution 
does not consider the possible influence of one or more 
other planets able to pump the eccentricity of WASP-6b 
(Mardling 2007) , but our RV data do not reveal the pres- 
ence of another planet so it seems reasonable at this stage 



to assume that the orbital evolution of WASP-6b was not 
dominated by planet-planet interactions. The model as- 
sumes also a constant radius for the planet during the 
whole tidal evolution, which is not very likely (Liu et al. 
2008). Furthermore, Fig. 6 shows that it does not conserve 
Ltot for e >^ 0.3 and during the final runaway merging 
of the planet with the star. Considering as valid only the 
part of the tidal evolution for which L tot is conserved at 
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Time (Gyr) Time (Gyr) 




Time (Gyr) Time (Gyr) 



Fig. 6. Tidal evolution for WASP-6b computed using the 
method described in Jackson et al. (2008b) . Top left: evo- 
lution of the semi-major axis. The dashed line shows the 
Roche limit of the system. Top right: evolution of the ec- 
centricity. Middle left: evolution of the tidal heating rate. 
Middle right: evolution of the total angular momentum. 
Bottom: evolution of the orbital (left) and star rotational 
(right) period. For each parameter, the solid line shows 
the evolution computed with the best-fitting present ec- 
centricity and semi-major axis while the dotted lines as- 
sume the maximum and minimum tidal heating consistent 
with their l-cr error bars. 

the 1-% level, we can nevertheless conclude from Fig. 6 
that WASP-6b experienced 0.6 - 1.2 Gyr ago a large tidal 
heating rate of 5 -10 x 10 19 W. Such a large heating rate in 
the past should have modified drastically the thermal his- 
tory of the planet and could have contributed significantly 
to the measured inflated radius. 

With a stellar irradiation ~ 4.7 10 8 erg s" 1 cm~ 2 , 
WASP-6b belongs to the theoretical pL planetary class 
proposed by Fortney et al. (2008; see also Burrows et al. 
2008). Under this theory, Ti and V-bearing compounds 
should mostly be condensed in the planetary atmosphere 
and secondary eclipse measurements at different wave- 
lengths should not reveal any stratospheric thermal in- 
version. Such secondary eclipses observations would not 
only constraint atmospheric models of giant close-in plan- 
ets, they would also constraint the eccentricity of the orbit 
and thus the tidal thermal history of the planet. 

The value that we determine for the sky-projected an- 
gle between the stellar spin and the planetary orbital axis 
is compatible with zero (j3 — lltig deg). This good align- 
ment was observed for ten other close-in giant planets (see 
Winn 2008 and references therein), while a possible mis- 



alignment was observed only for the planet XO-3 (Hebrard 
et al. 2008). Together, these results favor migration via 
tidal interactions with a protoplanetary disk (Lin et al. 
1996) as the dominant mechanism of planetary migration, 
because it should preserve spin-orbit alignement (Ward & 
Hahn 1994) contrary to migration via planet-planet scat- 
tering (Rasio & Ford 1996) or Kozai cycles (Fabrycky & 
Tremainc 2007). 

Acknowledgements. We thank the ESO La Silla staff for their 
support during the CORALIE and HARPS observations. The 
HARPS consortium is gratefully acknowledged for contributing 
to some of the observations presented here. 
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